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IMPINGEMENT  OF  WATER  DROPLETS  ON  A  RECTANGULAR  HALF  BODY  IN  A  TWO- 
DIMENSIONAL  INCOMPRESSIBLE  FLOW  FIELD 
By  William  Lewis  and  Rinaldo  J.  Brun 


SUMMARY 

Trajectories  of  water  droplets  moving  in  the  ideal  two-dimensional 
flow  field  ahead  of  a  hody  of  rectangular  cross  section  and  infinite  ex¬ 
tent  in  the  downstream  direction  have  been  calculated  by  means  of  a  dif¬ 
ferential  analyzer.  Data  on  collection  efficiency  and  distribution  of 
water  impingement  are  presented. 


INTRODUCTION 

As  part  of  a  comprehensive  research  program  dealing  with  the  prob¬ 
lem  of  icing  protection  for  aircraft,  an  investigation  of  the  impingement 
of  cloud  droplets  on  airfoils  and  other  aircraft  components  has  been  un¬ 
dertaken  by  the  NACA  Lewis  laboratory.  In  addition  to  specific  aircraft 
components  (refs.  1  to  3),  certain  general  aerodynamic  bodies  such  as 
cylinders  (ref.  4)  and  ellipsoids  of  revolution  (ref.  5)  have  been  studied 
because  various  aircraft  components  may  be  approximately  represented  by 
these  shapes.  Among  the  general  body  shapes  of  interest  is  that  of  a 
rectangular  body  having  a  flat  surface  facing  the  airstream.  Such  a  body 
may  be  used  as  an  approximation  to  the  shape  of  the  sensing  elements  of 
certain  instruments  used  to  measure  icing-cloud  properties. 

Cloud-droplet  size  and  size  distribution  are  sometimes  measured  by 
collecting  samples  of  cloud  droplets  on  small  transparent  slides  coated 
with  oil  or  some  water-sensitive  substance  (refs.  6  to  8) .  The  slides 
are  then  photographed  to  obtain  a  record  of  the  droplet  size  distribution. 
Various  types  of  sampling  devices  have  been  used  in  which  the  collecting 
element  is  usually  a  flat  surface  perpendicular  to  the  airstream.  The 
evaluation  of  the  cloud-droplet  size  distribution  from  samples  obtained 
in  this  way  requires  a  knowledge  of  the  relation  between  droplet  diameter 
and  collection  efficiency.  Droplet -impingement  data  applicable  to  a  body 
of  rectangular  cross  section  would  also  be  useful  in  estimating  the  col¬ 
lection  efficiency  of  the  central  portion  of  the  rotating-disk  icing-rate 
meter  (refs.  7  to  9)  and  the  rectangular  bars  used  to  measure  the  spray 
distribution  in  an  icing  wind  tunnel  (ref.  10) . 
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Trajectory  calculations  based  on  ideal  two-dimensional  flow  about  a 
ribbon  (ref*  ll)  ,  which  is  essentially  a  rectangle  of  zero  chord,  have 
been  used  as  an  approximation  to  obtain  the  impingement  on  rectangular 
bodies  of  finite  chord.  However,  the  ideal  flow  field  ahead  of  a  ribbon 
is  not  a  very  satisfactory  representation  of  the  actual  flow  field  in 
front  of  a  rectangular  body  because  of  the  dependence  of  the  flow  ahead 
of  the  body  on  conditions  existing  behind  the  body. 

The  ideal  flow  field  about  a  ribbon  is  symmetrical  fore  and  aft,  as 
shown  in  figure  l(a) ,  with  the  streamlines  closing  in  directly  behind  the 
body.  In  the  case  of  real  flow  about  a  rectangular  body,  on  the  other 
hand,  the  streamlines  cannot  be  calculated  exactly;  but  general  knowledge 
of  the  flow  about  solid  boundaries  suggests  that  the  flow  field  is  some¬ 
thing  like  that  shown  schematically  in  figure  l(b) .  Flow  separation  oc¬ 
curs  at  the  edge,  and  a  wake  is  formed  extending  downstream.  The  effect 
of  variations  in  the  ratio  of  chord  to  thickness  cannot  be  estimated  ac¬ 
curately,  but  is  probably  small,  since  the  wake  is  roughly  equivalent  to 
an  indefinite  downstream  extension  of  the  body. 

Potential  flow  exists  outside  the  limits  of  the  boundary  layer  and 
wake  (shown  approximately  by  the  dashed  line).  Figure  1(c)  shows  the 
ideal-flow  streamlines  around  a  rectangular  half  body  (a  body  of  infinite 
extent  downstream)  .  This  flow  field  resembles  the  estimated  real  flow  in 
figure  1(b)  much  more  closely  than  does  the  ideal  flow  about  a  ribbon 
(fig.  1(a)).  The  effect  of  the  downstream  extent  of  the  body  on  the  flow 
field  ahead  is  shown  in  figure  l(d),  which  is  a  superposition  of  portions 
of  the  fields  shown  in  figures  l(a)  and  (c) .  The  comparison  shows  that 
the  presence  of  the  afterbody  exerts  a  considerable  influence  on  the  flow 
field  ahead  of  the  body. 

The  ideal  flow  field  of  the  rectangular  half  body  represents  a  closer 
approximation  to  real  flow  than  does  the  ideal  flow  about  a  ribbon;  how¬ 
ever,  it  too  is  subject  to  the  disadvantage  that  the  velocity  becomes  in¬ 
finite  at  the  corner.  In  spite  of  this  disadvantage,  the  ideal  flow  field 
in  front  of  the  rectangular  half  body  was  used,  in  the  calculation  of  tra¬ 
jectories,  because  it  was  considered  to  be  sufficiently  realistic  to  yield 
approximate  results  of  practical  value  and  because  the  velocity  components 
could  be  calculated  without  too  much  difficulty. 


SYMBOLS 

The  following  symbols  are  used  in  this  report: 
A  area  per  unit  span,  sq  ft/ft 

a  droplet  radius,  ft 

B,C  empirical  constants,  dimensionless 
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ux>uy 


vx'vy 


coefficient  of  drag,  dimensionless 

collection  efficiency,  dimensionless 
functional  relation  between  y q  and  yg 

flux  of  liquid  water  per  unit  span,  slugs/(sec) (ft) 
inertia  parameter,  dimensionless 
half  width  of  rectangle,  ft 

♦ 

impingement  distribution  index,  dimensionless 
arbitrary  numbers 

local  Reynolds  number  with  respect  to  droplet,  dimensionless 
free-stream  Reynolds  number  with  respect  to  droplet,  dimensionless 

time,  sec 

free-stream  velocity,  ft/sec 

local  air  velocity  components,  ratio  to  free-stream  velocity, 
dimensionless 

local  droplet  velocity  components,  ratio  to  free-stream  velocity, 
dimensionless 

total  rate  of  water  impingement  per  unit  span,  slugs/ (sec)  (ft) 
local  rate  of  water  impingement  per  unit  area,  slugs/ ( sec )  (sq  ft) 

free-stream  liquid- water  content,  slugs/cu  ft 
rectangular  coordinates,  ratio  to  half  width  L,  dimensionless 
distance  from  centerline  (Y  =  Ly) ,  ft 
trajectory  starting  ordinate,  dimensionless 

trajectory  ordinate  at  point  of  impingement,  dimensionless 


local  impingement  efficiency,  dimensionless 
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Po  local  impingement  efficiency  at  centerline 

P]_  local  impingement  efficiency  at  edge 

H  viscosity  of  air,  slugs/ (ft) (sec) 

Pa  density  of  air,  slugs/cu  ft 

density  of  water,  slugs/cu  ft 

t  dimensionless  time  function,  t  =  — 

*  L 

®  Re^/K,  dimensionless 


ANALYSIS 

In  order  to  find  the  rate  and  distribution  of  droplet  impingement 
on  the  surface  of  a  body,  it  is  necessary  to  determine  the  cloud-droplet 
trajectories  with  respect  to  the  body.  The  method  used  in  this  case  to 
calculate  the  trajectories  of  cloud  droplets  is  described  fully  in  ref¬ 
erence  4.  Assumptions  that  are  necessary  to  the  solution  of  the  problem 
are:  (l)  At  a  large  distance  ahead  of  the  body  the  droplets  are  at  rest 

with  respect  to  the  air;  (2)  the  only  external  force  acting  on  the  drop¬ 
lets  is  the  drag  force  due  to  the  relative  velocity  of  the  air  with  re¬ 
spect  to  the  droplets;  and  (3)  the  droplets  are  always  spherical  and  do 
not  change  in  size. 


Differential  Equations  of  Droplet  Motion 

The  differential  equations  that  describe  the  motion  of  the  droplets 
were  obtained  by  equating  the  drag  force  with  the  rate  of  change  of  mo¬ 
mentum  of  the  droplet  (ref.  4).  The  equations  are  expressed  in  dimension¬ 
less  form  in  order  to  maintain  the  number  of  calculations  at  a  minimum 
and  to  simplify  the  presentation  of  the  results.  They  apply  to  the  mo¬ 
tion  of  droplets  in  a.  plane  perpendicular  to  the  edges  of  the  rectangular 
half  body  in  the  system  of  rectangular  coordinates  shown  in  figure  2. 

The  origin  is  located  at  the  center  of  the  front  face  of  the  rectangle, 
and  the  unit  of  distance  is  the  half  width  L  of  the  rectangle.  The 
dimensionless  coordinates  x  and  y  are  ratios  to  L.  The  dimension¬ 
less  air  velocity  components  u^  and  u  and  droplet  velocity  components 

y 

vx  and  vy  are  ratios  to  the  free-stream  velocity  U,  which  is  the  unit 

of  velocity.  Time  is  expressed  in  terms  of  the  dimensionless  quantity 
t  =  tU/L.  The  unit  of  time  is  the  time  required  to  go  a  distance  L  at 
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a  speed  U.  At  a  large  distance  ahead  of  the  body,  a  uniform  air  flow 
bearing  cloud  droplets  is  assumed  to  be  approaching  the  rectangle  from 
the  negative  x  direction  moving  parallel  to  the  x  axis  at  a  speed  U 
(fig.  2). 

The  equations  of  motion  of  the  droplets  are 


dvx  _  CjjRe  x 

dx  24  K 

i  (u  -  v ) 

dx  24  K  '  y  y' 


where  Re  is  the  local  Reynolds  number  with  respect  to  the  droplet  diame¬ 
ter  2a  and  the  local  relative  velocity  between  the  air  and  the  droplet, 
thus 

Re  =  -  vx)2  +  (Uy  -  vy)2  (2) 

The  coefficient  in  equation  (2)  is  called  the  free -stream  Reynolds  number 
ReQ;  hence, 

2ap  U 

_  Of  /  *7  \ 


The  dimensionless  number  K,  the  inertia  parameter,  is  defined  as  follows 

K  =  (4) 

The  coefficient  of  drag  CD  may  be  obtained  from  experimental  data  as  a 
function  of  the  local  droplet  Reynolds  number  Re. 

An  examination  of  equations  (l)  and  (2)  shows  that  the  characteris¬ 
tics  of  the  trajectories  depend  only  on  the  values  of  the  two  dimension¬ 
less  parameters  K  and  Reg*  Thus,  a  unique  set  of  droplet  trajectories 

exists  for  every  combination  of  K  and  Reg. 


Method  of  Solution 

The  differential  equations  of  motion  (l)  are  difficult  to  solve  be¬ 
cause  values  of  the  velocity  components  and  the  factor  containing  the 
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coefficient  of  drag  depend,  on  the  position  and  velocity  of  the  droplet  at 
each  instant  and>  therefore^  are  not  known  until  the  trajectory  is  traced. 
The  values  of  these  quantities  must  "be  fed  into  the  equations  as  a  tra¬ 
jectory  is  developed.  This  was  accomplished  by  using  a  mechanical  dif¬ 
ferential  analyzer  constructed  at  the  NACA  Lewis  laboratory  for  this  pur¬ 
pose  (ref.  4)  .  The  results  were  obtained  in  the  form  of  plots  of  droplet 
trajectories  in  the  coordinates  of  figure  2. 

The  following  information  was  required  for  use  with  the  differential 
analyzer  in  the  trajectory  calculation: 

(1)  Coefficient  of  drag  data:  Values  of  the  factor  cD  Re/24  as  a 
function  of  Re  were  obtained  from  tables  in  reference  11. 

(2)  Air  velocity  components:  The  air  velocity  components  u  and 

X 

Uy  were  determined,  as  functions  of  x  and  y  by  means  of  an  analytical 

solution  of  the  two-dimensional  potential  flow  field  ahead  of  a  rectangu¬ 
lar  half  "body.  Because  of  symmetry  with  respect  to  the  centerline  (see 
^■8*  }  one  half  of  this  field  is  identical  to  the  flow  field  around  a 

step  discontinuity  in  the  floor  of  a  channel  of  infinite  depth,  a  configu¬ 
ration  discussed  in  reference  12.  The  equations  used  to  calculate  the 
velocity  components  were  obtained  by  means  of  the  Schwarz-Christoffel 
transformation  as  described  in  reference  12.  These  equations,  expressed 
in  parametric  form,  are  as  follows: 


Any  pair  of  values  of  the  arbitrary  parameters  p  and  q  deter¬ 
mines  the  coordinates  x  and  y  of  a  point  in  the  flow  field  and  the 
components  u  and  u  of  the  velocity  at  that  point.  Values  of  u  . 

Uy,  x,  and  y  were  calculated  from  equation  (5)  for  a  number  of  points 

in  the  flow  field  ahead  of  the  rectangle,  and  curves  giving  u^.  and  uy 

as  functions  of  x  and  y  were  obtained  by  cross-plotting.  The  calcu¬ 
lated  values  of  coordinates  and  air  velocity  components  are  presented  in 
table  I . 
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(3)  Starting  conditions:  The  integrations  were  started  at  x  =  -60 
where  the  air  velocity  components  are  u^.  =  0.9945  and  Uy  =  O.OOOly 

for  all  values  of  y  between  0  and  1.  The  initial  droplet  velocity 
components  were  assumed  to  he  vx  =  0.9945  and  vy  =  0.  The  calcula¬ 
tions  were  performed  in  three  steps:  x  =  -60  to  -10;  x  =  -10  to  -1; 
and  x  =  -1  to  0.  Different  scale  factors  were  used  on  the  machine  for 
each  step,  and  the  results  of  each  of  the  first  two  steps  were  used  to 
establish  starting  conditions  for  the  succeeding  step. 


Determination  of  Droplet  Impingement 

The  differential  analyzer  traces  the  trajectories  of  the  droplets 
from  a  starting  point  defined  by  the  ordinate  yQ  at  a  large  distance 

ahead  of  the  body  to  the  final  point  of  impingement  having  the  ordinate 

y  at  the  surface  of  the  rectangle.  The  trajectories  provide  pairs  of 
s 

values  of  yQ  and  yg  that  establish  the  relation  between  yQ  and  ys 
for  various  combinations  of  K  and  RSq.  The  relation  is  expressed  sym¬ 
bolically  as  follows: 

y0  -  f (y»)  <6> 


and  the  symbols  f(o)  and  f (l)  are  used  to  denote  the  values  assumed  by 

y„  when  y=  is  0  and  1,  respectively, 
u  s 

If  a  rectangular  body  of  half  width  L  moves  through  a  cloud  of 
uniform  liquid- water  content  and  uniform  droplet  size,  the  flux  of  liquid 
water  per  unit  span  between  any  two  trajectories  with  starting  ordinates 
LyQ  ^  and  Ly0  2  (LyQ  has  the  dimensions  of  length)  is 

G  =  UwL(y0^2  -  yo,l) 


Since  the  trajectories  define  the  path  of  the  flow  of  liquid  water,  the 
same  flux  arrives  at  the  rectangle  surface  between  the  trajectory  terminal 
ordinates  Ly  n  and  Ly  _  corresponding  to  the  starting  ordinates 

S  •  X  S  a  Cx 


Ly0,i  and  LyQ^. 
rectangle  between 


The  surface  area  per  unit  span  of  the  strip  of  the 
the  two  terminal  ordinates  Lyg  1  and  Lyg  2  is 


A  =  L<ys,2  “  y! 


The  liquid-water  flux  G  is  distributed  over  the  area  A;  hence,  the 
average  rate  of  water  impingement  per  unit  area  on  this  strip  is 


a 
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e  ""^o.z  -  r0.i> 

A  ~  Liys,2  -  yB,i) 

If  ys,l  is  allowed  to  approach  ys^2  as  a  limit,  the  limit  approached 
by  G/A  is  the  local  rate  of  water  impingement  W  at  the  distance 

y 

Y  =  ygL  from  the  centerline. 


G  dy0 

W„  -  lim  f  =  Uw  ~ 

y  A-»0  A 


The  local  impingement  efficiency  (3  is  defined  as  follows 


3 


Hence , 


3  = 


dY0 


(y) 


(9) 


The  total  rate  of  water  interception  per  unit  span  over  the  entire 
face  of  the  rectangle  is  twice  the  integral  of  the  local  impingement  rate 
from  the  stagnation  point  (Y  =  0)  to  the  outside  edge  (Y  =  L);  thus. 


W  =  2LUw[f(l)  -  f(0)J 


where  the  symbol  f  refers  to  the  function  of  equation  (6) . 

In  a  symmetrical  flow  field,  a  droplet  approaching  the  stagnation 
point  along  the  central  streamline  does  not  change  direction;  thus,  the 
value  of  Yq  corresponding  to  ys  =  0  is  0;  therefore. 


W  =  2LUwf  (l) 

The  collection  efficiency  E  is  defined  as  follows: 


(10) 


2LUw 


(ID 


3804 


3804  2-10 


NACA  TN  3658 


9 


Hence , 

E  =  f (1)  (12) 

The  collection  efficiency  is  therefore  equal  to  the  value  of  y q  when 
yg  =  1,  which  is  the  nondimensional  starting  ordinate  of  the  trajectory 

that  strikes  the  outside  edge  of  the  rectangle  in  the  coordinate  system 
of  figure  2. 


RESULTS 

A  total  of  74  pairs;  of  values  of  yQ  and  ys  were  calculated  with 
the  differential  analyzer  using  23  comhinations  of  ReQ  and  K.  It  was 
found  that  the  relation  of  equation  (6)  between  yQ  and  ys  for  any 
particular  combination  of  K  and  ReQ  could  be  represented  approximate¬ 
ly  by  an  empirical  equation  of  the  form 

y0  =  ^s  +  cys  •  (13) 

where  B  and  C  are  positive  and  less  than  1  and  n  is  positive  and 
greater  than  1.  The  collection  efficiency  E  is  the  value  of  yQ  when 

ys  =  1;  hence , 

E  =  B  +  C  (14) 

The  local  impingement  efficiency  0  may  be  obtained  by  differenti¬ 
ating  equation  (13)  as  follows: 

p  =  =  B  +  nCy^  (15) 

"  § 

The  local  impingement  efficiency  at  the  stagnation  point  pQ  is  obtained 
by  setting  yg  =  0  in  equation  (15) j  thus. 


and  from  equation  (14), 


B  =  P, 


C  =  E  -  P, 


Substituting  in  equation  (13)  gives 


yo  -  0oys  +  (b  -  0o)ys 


(16) 
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and  from  equation  (15)  there  is  obtained 

P  =  0O  +  n(E  -  P0)y(n“l)  (17) 

Equation  (17)  shows  that  the  local  impingement  efficiency  has  a  minimum 
value  (3q  at  the  stagnation  point  and  increases  to  a  maximum  given  by 

Pl  =  P0  +  n(E  -  P0)  (18) 


at  the  outside  edges  of  the  rectangle. 

The  impingement  pattern  corresponding  to  any  given  pair  of  values 
of  K  and  Re^  is  determined  by  three  parameters:  (l)  the  collection 

efficiency  E,  (2)  the  local  impingement  efficiency  at  the  stagnation 
point  P0,  and  (3)  the  distribution  index  n,  which  determines  the  form 

of  the  distribution  of  impingement. 

An  analysis  of  the  trajectory  data  has  shown  that  these  three  param¬ 
eters,  which  are  functions  of  K  and  ReQ,  may  be  represented  approxi¬ 
mately  by  the  following  empirical  equations: 


.1 

n  =  1  +  e  K  +  5eQ-+°52  t1*1  +  (loee0*72  K)2]"2  (19) 


l0SeP°  =  (kTTo  +  0,2  +  °'0238  Reo'538)(] 


1.16  - 


115  +  Re 


87 _ \ 

q  +  0 • 004  R6q J 


(K  -  0.28) 

151 


-0.6 


(20) 


-logeE  - 


+  150 


+  0.267  +  0.225  Re 


02 


180 


250  +  Re  q 


) 


(K  -  0.15) 


-0.74 


(21) 


Equations  (17) ,  (19),  (20),  and  (21)  provide  a  means  of  determining 
approximately  the  collection  efficiency  and  the  distribution  of  impinge¬ 
ment  for  any  combination  of  ReQ  and  K.  The  trajectory  calculations 

used  in  deriving  the  equations  adequately  covered  a  range  of  K  from 
l/2  to  10  and  of  ReQ  from  16  to  256.  Outside  this  range  the  equa¬ 
tions  represent  extrapolations. 
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Values  of  n,  0Q,  and  E  calculated  from  the  empirical  equations 

are  presented  in  figures  3,  4,  and  5,  respectively.  The  74  values  of 
yQ  obtained  with  the  differential  analyzer  and  the  corresponding  values 

calculated  from  the  equations  are  presented  for  comparison  in  table  II. 

The  standard  deviation  of  the  differences  between  corresponding  values 
is  0.0051;  therefore ,  the  probable  error  involved  in  using  the  equations 
to  determine  yQ  is  ±0.0034  based  on  the  assumption  that  the  results  from 

the  trajectories  are  absolutely  correct.  The  dimensionless  parameters 
ReQ  and  K  are  defined  in  equations  (3)  and  (4) .  Equations  and  graphs 

for  use  in  determining  values  of  K  and  ReQ  in  terms  of  air  speed,  al¬ 
titude,  body  size,  and  droplet  diameter  in  practical  units  are  presented 
in  reference  13. 


DISCUSSION 

The  use  of  empirical  equations  to  represent  the  impingement  patterns 
is  essentially  a  means  of  interpolating  and  extrapolating  the  trajectory 
data.  The  procedure  is  analogous  to  that  of  drawing  a  family  of  curves 
and  interpolating  between  them.  Because  of  random  errors  in  the  trajec¬ 
tory  data,  a  certain  amount  of  smoothing  is  necessary  in  using  either 
procedure.  The  use  of  empirical  equations  has  the  advantage  of  accomplish¬ 
ing  the  smoothing  operation  in  a  consistent  manner  over  the  entire  range 
of  K  and  Re„,  thus  assuring  that  the  results  are  as  accurate  for  inter¬ 
mediate  values  of  these  parameters  as  for  the  values  for  which  trajecto¬ 
ries  were  calculated.  The  equations  also  provide  a  more  reliable  basis 
for  extrapolation  to  values  of  K  and  ReQ  beyond  the  range  of  the 
trajectory  calculations. 

The  average  accuracy  of  the  equations  is  probably  a  little  better 
than  is  Indicated  by  the  probable  error  of  ±0.0034  in  yQ  based  on  the 

departures  in  table  II,  because  these  deviations  are  due  partly  to  errors 
in  the  trajectory  calculations.  The  same  degree  of  accuracy  (probable 
error  about  0.003)  also  applies  to  values  of  E  and  PQ  from  equations 

(20)  and  (2l).  Values  of  3  at  intermediate  points  and  values  of  0-l 

at  the  corner  of  the  rectangle,  calculated  from  equations  (17)  and  (18), 
may  be  subject  to  larger  errors  because  of  the  effect  of  errors  in  the 
distribution  index  n. 

The  validity  of  the  representation  of  the  impingement  distribution 
patterns  given  by  equation  (17),  using  values  of  n,  30,  and  E  from 

equations  (19),  (20),  and  (2l),  respectively,  can  be  checked  by  a  com¬ 
parison  with  average  values  of  3  for  the  intervals  of  ys  between 

adjacent  trajectories.  Such  values  are  obtained  directly  from  the  tra¬ 
jectory  data  by  means  of  the  relation 
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The  comparison  is  shown  in  figure  6  in  which  the  smooth  curves  were 
calculated  from  equation  (17),  and  the  block  diagrams  represent  averages 
from  equation  (22)  based  on  individual  intervals  between  adjacent  cal¬ 
culated  trajectories.  An  examination  of  figure  6  reveals  the  general 
pattern  of  the  impingement  distribution  as  a  function  of  K  and  ReQ. 

At  high  values  of  K,  the  impingement  efficiency  approaches  1  across  the 
entire  face  of  the  rectangle,  but  more  slowly  for  high  values  of  ReQ. 

At  low  values  of  K,  the  impingement  efficiency  approaches  0,  but  the 
number  of  trajectories  calculated  (3  at  K  =  l/3)  is  too  small  to  deter¬ 
mine  accurately  the  shape  of  the  impingement  profile.  For  all  values  of 
K,  the  impingement  efficiency  decreases,  and  the  distribution  becomes 
more  uniform  with  increasing  ReQ.  The  greatest  curvature  of  the  impinge¬ 
ment  profile,  indicating  a  sharp  maximum  impingement  at  the  edge,  occurs 
at  low  values  of  ReQ  and  values  of  K  between  1  and  2.  The  curves  cal- 

culated  from  equation  (17)  agree  satisfactorily  with  the  values  of  "p  ob¬ 
tained  from  adjacent  trajectories  over  the  entire  range  of  K  and  ReQ 

for  which  trajectories  were  determined.  Equation  (17),  therefore,  pro¬ 
vides  a  consistent  representation  of  the  variations  of  the  impingement 
profile  with  K  and  ReQ. 

It  is  of  interest  to  compare  the  droplet-impingement  data  for  a 
rectangle  with  similar  data  for  a  ribbon  in  reference  11.  The  results 
are  presented  in  reference  11  in  the  form  of  curves  showing  E  and  p 
as  functions  of  K  and  ®  (d»  =  Reg/ft) .  Values  of  E  and  pQ  read 
from  these  curves  are  compared  in  table  III  with  values  calculated  from 
equations  (16)  and  (17)  for  corresponding  values  of  K  and  ReQ.  For 

most  values  of  K  and  $  the  rectangle  has  a  lower  collection  effi¬ 
ciency  than  the  ribbon;  this  effect  is  due  to  the  influence  of  the  after¬ 
body  in  reducing  the  rate  of  deceleration  of  the  air  in  the  region  iust 
ahead  of  the  body. 

It  was  pointed  out  in  the  discussion  of  flow  fields  in  the  INTRODUC¬ 
TION  that  the  flow  field  of  the  rectangular  half  body  differs  from  that 
of  the  ribbon  by  the  addition  of  an  afterbody,  but  that  both  differ  from 
the  real  flow  about  a  rectangular  body  because  of  the  absence  of  flow 
separation  and  because  of  the  requirement  of  infinite  velocity  at  the 
edges.  An  examination  of  the  comparative  streamlines  in  figure  1(d)  shows 
a  great  difference  in  velocity  in  the  immediate  vicinity  of  the  edge. 
(Velocity  is  inversely  proportional  to  streamline  spacing.)  In  spite  of 
the  difference  in  the  flow  fields,  however,  the  differences  in  collection 
®  ficiency  of  the  ribbon  and  rectangle  are  only  moderate.  The  maximum 
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difference  in  table  III  is  0.099  at  K  =  2  and  $  =  10,000,  and  the 
average  difference  for  all  values  in  table  III  is  0.026.  These  facts 
suggest  that, the  results  obtained  using  the  ideal  flow  ahead  of  the  rec¬ 
tangular  half  body  are  likely  to  be  a  fairly  good  approximation  to  the 
actual  impingement  on  a  rectangular  body  in  real  flow.  The  largest  dif¬ 
ference  between  the  ideal  flow  and  the  real  flow  is  that  the  ideal  flow 
permits  a  high  local  air  velocity  in  the  neighborhood  of  the  edge;  how¬ 
ever,  these  differences  in  local  air  velocity  in  the  neighborhood  of  the 
edge  have  a  relatively  small  effect  on  the  trajectories. 


CONCLUDING  REMARKS 

The  rate  and  distribution  of  water  impingement  on  a  semi -infinite 
rectangular  body  moving  in  a  cloud  of  uniform  liquid-water  content  and 
droplet  size  have  been  determined  by  means  of  droplet  trajectory  calcula¬ 
tions.  Empirical  equations  have  been  found  providing  an  approximate  rep¬ 
resentation  of  the  impingement  rate  and  distribution  in  terms  of  the 
inertia  parameter  K  and  the  free-stream  droplet  Reynolds  number  ReQ. 

Equation  (17)  gives  the  local  impingement  efficiency  (3  as  a  function  of 
distance  from  the  centerline  of  the  rectangle  in  terms  of  the  following 
three  parameters:  (l)  the  collection  efficiency  E,  (2)  the  impingement 
efficiency  at  the  centerline  p  ,  and  (3)  the  impingement  distribution 

index  n. 


These  parameters  may  be  determined  from  K  and  ReQ  by  means  of 
empirical  equations  (19),  (20),  and  (2l) . 


Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  December  14,  1955 
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TABLE  I.  -  VELOCITY  FIELD  AHEAD  OF  RECTANGULAR  HALF  BODY 

[The  table  gives  values  of  the  dimensionless  space 
coordinates  x  and  y  and  corresponding  values 
of  dimensionless  velocity  components  u^.  and 

Uy  (see  fig.  2)  calculated  from  equation  (5)]| 


X 

y 

“x 

u 

y 

X 

y 

”x 

uy 

-0.09567 

0 

0.07485 

0 

0 

0.12551 

0 

0.09918 

-.19242 

.14889 

-.09473 

.12622 

.07571 

.09855 

-.29136 

. 22128 

-.19053 

.12835 

.15036 

.09699 

-.39364 

.29132 

- . 28847 

.13195 

. 22335 

.09429 

-.50049 

.35837 

-.38969 

.13712 

. 29394 

.09073 

-.61323 

.42191 

-.49538 

.14395 

.3615 

.08637 

-.73335 

.48154 

- .  60686 

.15260 

.4253 

.08142 

-.86246 

.53707 

-.72556 

.16328 

.4852 

.07603 

-1.0518 

.60437 

?* .  85310 

.17622 

.5408 

.07040 

-1.2655 

. 66404 

-1.0400 

.19750 

.6081 

.06275 

-1.6451 

.73977 

-1.2507 

. 22419 

.6677 

.05523 

-2.1190 

. 80050 

-1.6248 

.27674 

.7432 

.04469 

-2.9586 

.86173 

-2.0914 

.34869 

.8033 

.03542 

-4.1437 

.90515 

-2.9171 

.48638 

.8639 

.02536 

-6.3798 

.94138 

-4.0816 

.69271 

.9067 

.01779 

-9.9599 

.96403 

-6.2771 

1.0996 

.9424 

.01119 

-11.96 

.97045 

-9.7905 

1.7711 

.9647 

.00695 

-14.40 

.97580 

0 

. 24681 

0 

.19892 

-17.30 

.98010 

-.09202 

. 24818 

.07782 

.19774 

-20.9 

.98368 

-. 18505 

.25233 

.15466 

.19434 

-25.2 

.98660 

- . 28009 

. 25935 

. 22960 

.18880 

-30.5 

.98903 

-.37821 

.26940 

.30184 

.18144 

-37.0 

.99100 

-.48056 

.28271 

.3707 

.17251 

-44.8 

.99264 

-.58836 

.29957 

.4356 

.16237 

-54.5 

.99397 

-.70299 

.32036 

.4960 

.15138 

-66.0 

.99505 

-.82589 

.34556 

.55204 

.13995 

-98 

.99668 

-1.00567 

.38700 

.61937 

.12446 

-145 

.99777 

-1.2079 

.43897 

.67846 

.10929 

-216 

.99851 

i 

-1.5659 

.54132 

.75277 

.08814 

-320 

.  1 

1 

.99900  T 

-2.0110 

.68141 

.81161 

.06968 

CM  o 
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TABLE  I .  -  Continued.  VELOCITY  FIELD  AHEAD  OF  RECTANGULAR  HALF  BODY 

[The  table  gives  values  of  the  dimensionless  space 
coordinates  x  and  y  and  corresponding  values 
of  dimensionless  velocity  components  u^.  and 

Uy  (see  ■£>ig*  2)  calculated  from  equation  (5) .] 


x 


-2.7964 

-3.9010 

0 

-.08763 

-.17617 

-.26652 

-.35965 

-.45658 

-.55842 

-.66641 

- . 78186 
-.95012 
-1.1386 
-1.4706 
-1.8804 

6012 

-.08163 
- . 16403 
- . 24798 

-.33428 

-.42382 

-.51753 

-.61646 

-.70243 

-.87424 

-1.0440 

-1.3404 

0 

-.07431 


y 


0.94955 

1.3514 

.36350 

.36549 

.37150 

.38167 

.39623 

.41550 

.43993 

.47005 

.50656 

.56659 

.64187 

.79014 

.99310 

1.3815 

.47508 

.47762 

.48529 

.49826 

.51683 
.54141 
.57258 
. 61100 
.65756 

.73412 

.83015 

1.0193 

.57788 

.58087 


0.87027 

.91143 

0 

.08318 

.16220 

.24046 

.31554 

.38666 

.45316 

.51469 

.57037 

.63835 

.69669 

.76902 

.82538 

.88073 

0 

.08763 
.17377 
. 25708 

.33641 

.41084 

.47965 

.54257 

.59954 

.66631 

.72327 

.79232 

0 

.09618 


u 

y 


0.04975 

.03482 

.30257 

.30063 

.29528 

.28647 
.27476 
. 26064 
.24473 
. 22756 

. 20947 
.18583 
.16258 
.13044 
.10266 

.07293 

.41422 

.41151 

.40350 

.39066 

.37364 
.35325 
.33043 
.30599 
. 28086 

. 24740 
.21528 
.17142 
.53542 
.53158 


x 


-0.14923 

-.22535 

-.30333 

-.38384 

-.46762 

_- .  55549 
-.64834 
-.78167 
-.92848 
0 

-.06609 
-.13259 
-.19992 
- . 26854 
-.33890 

-.41152 

-.48696 

-.56585 

-.67760 

0 

-.05716 

-.11454 

-.17235 

-.23081 

-.29017 

-.35069 

-.41266 

-.47639 

0 

-.04775 


y 


0.58990 

.60516 

.62702 

.65597 

.69261 

.73784 

.79265 

.88280 

.99580 

.66908 

.67239 

.68241 

.69936 

.72361 

.75573 

.79644 
.84664 
.90750 
1.0075 
[  .74914 

.75266 

.76331 

.78129 

.80704 

.84113 

.88433 

.93758 

1.0022 

.81831 

.82190 


”x 


0.19037 

.25295 

.36594 

.44476 

.51652 

.58097 

.63827 

.70392 

.75851 

0 

.10801 
. 21326 
.31323 
.40600 
.49028 

.56537 

.63121 

.68822 

.75165 

0 

.12466 

.24522 

.35810 

.46069 

.55148 

.62994 

.69645 

.75206 

0 

.14888 


°y 


0.52027 

.47837 

.45015 

.41876 

.38561 

.35190 

.30766 

.26574 

.66820 

.66394 

.64699 

.62186 

.58914 

.55077 

-  .50883 
.46506 
.42123 
.36466 
.81949 

.81181 

.78960 

.75456 

.70953 

.65756 

.60170 

.54458 

.48848 

1.0000 

.98885 


* 


CI-3  3804 


MCA  TN  3658 


17 


TABLE  I.  -  Concluded.  VELOCITY  FIELD  AHEAD  OF  RECTANGULAR  HALF  BODY 


[The  table  gives  values  of  the  dimensionless  space 
coordinates  x  and  y  and  corresponding  values 
of  dimensionless  velocity  components  u^.  and 

Uy  (see  fig.  2)  calculated  from  equation  (5).] 


X 

y 

X 

y 

ux 

"y 

-0.09549 

0.83274 

0.29131 

0.95663 

-0.02118 

0.95495 

0.32951 

1.8215 

-.14324 

.85109 

.42190 

.90665 

-.04176 

.96398 

.62023 

1.6952 

-.19099 

.87735 

.53705 

.84356 

-.06113 

.97924 

.76417 

1.3668 

-.23873 

.91212 

.63516 

.77240 

-.07864 

1.0011 

1.0088 

1.3188 

-.28648 

.95617 

.71632 

.69781 

0 

.97508 

0 

2.4142 

-.33423 

1.0105 

.78179 

.62349 

-.01386 

.97762 

.49501 

2.3248 

0 

.87506 

0 

1.2203 

-.02696 

.98529 

.90034 

2.0906 

-  .03833 

.87858 

.18480 

1.2034 

-.03850 

.99826 

1.1755 

1.7863 

-.07645 

.88923 

.35876 

1.1552 

0 

.98942 

0 

3.3050 

-.11413 

.90721 

.51336 

1.0817 

-.00786 

.99141 

.85669 

3.1543 

-.15116 

.93296 

.64377 

.99151 

-.01482 

.99742 

1.4292 

2.6018 

-.18729 

.96705 

.74884 

.89289 

0 

.99681 

0 

5.0275 

-.22227 

1.0103 

.83013 

.79296 

-.00347 

.99818 

1.7229 

4.3791 

0 

.91908 

0 

1.4966 

0 

.99959 

0 

10.078 

-.02940 

.92239 

.23954 

1.4687 

-.00076 

1.0003 

4.8931 

6.3865 

-.05840 

.93241 

.45956 

1.3943 

-.08656 

.94936 

.64606 

1.2826 

-.11343 

.97361 

.79306 

1.1508 

-.13856 

1.0057 

.90171 

1.0130 

0 

.95196 

0 

1.8675 
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TABLE  II.  -  COMPARISON  OF  CALCULATED  STARTING  ORDIMTES 

[Starting  ordinates  yQ  of  trajectories  calculated  with  differential 
analyzer  and  corresponding  values  of  yQ  calculated  from  empiri¬ 
cal  equations  for  same  values  of  ys  are  compared.] 


yo  yo 

(Trajectory)  (Equations) 


0.676 

1.028 


0.414 

.551 

.807 

1.012 


0.238 

.472 

.703 

.926 


0.425 

.999 


0.333  0.414 


0.370 

.554 

.755 

.904 


0.228 

.461 


0.033 

.062 

0.075 

.149 

.224 

.265 

0.082 

.163 

.245 

.327 

.408 

.492 

.505 


0.259 

.346 

.523 

.702 


0.186 

.369 

.557 

.740 


0.384 

.916 


0.013 


0.037 

.054 

.088 

.109 


0.227 

.461 


0.037 

.070 

0.074 

.149 

.225 

.257 

0.082 

.166 

.250 

.327 

.405 

.495 

.503 


0.266 

.355 

.528 

.697 


0.186 

.371 

.557 

.745 

0.384 

.916 


0.012 


0.037 

.062 

.095 

.124 


Difference 


-0.001 

0 


+0.004 
+  .008 


-0.001 

0 

+.001 

-.008 


0 

+  .003 
+  .005 
0 

-.003 
+  .003 
-.002 


+0.007 
+  .009 
+  .005 
-.005 


0 

+  .002 

0 

+  .005 


-0.001 


0 

+  .008 
+  .007 
+  .015 
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TABLE  II.  -  Continued.  COMPARISON  OF  CALCULATED  STARTING  ORDINATES 

[Starting  ordinates  yQ  of  trajectories  calculated  with  differential 
analyzer  and  corresponding  values  of  Yq  calculated  from  empiri¬ 
cal  equations  for  same  values  of  ys  are  compared.] 


K 

yo 

yo 

Difference 

(Trajectory) 

(Equations) 

1 

0.127 

0.040 

0.037 

-0.003 

.490 

.158 

.146 

-.012 

.722 

.237 

.229 

-.008 

.974 

.359 

.352 

-.007 

1.035 

.400 

.390 

-.010 

2 

0.500 

0.251 

0.254 

+0.003 

.644 

.338 

.333 

-  .005 

.920 

.511 

.511 

0 

4 

0.261 

0.180 

0.178 

-0.002 

.519 

.361 

.361 

0 

.765 

.542 

.545 

+  .003 

.985 

.727 

.720 

-.007 

10 

0.443 

0.378 

0.378 

0 

.977 

.857 

.856 

-.001 

0.5 

0.545 

0.035 

0.034 

-0.001 

.820 

.053 

.060 

+  .007 

1 

0.184 

0.038 

0.037 

-0.001 

.351 

.077 

.072 

-.005 

.515 

.114 

.110 

-.004 

.638 

.153 

.143 

-.010 

.910 

.232 

.233 

+  .001 

2 

0.399 

0.164 

0.164 

0 

.576 

.248 

.244 

-.004 

.748 

.331 

.331 

0 

.931 

.440 

.436 

-  .  004 

4 

0.287 

0.176 

0.174 

-0.002 

.565 

.352 

.354 

+  .002 

.700 

.440 

.446 

+  .006 

.944 

.621 

.621 

0 
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TABLE  III.  -  COMPARISON  OF  COLLECTION  EFFICIENCY  E  AND  STAGNATION 


LINE  IMPINGEMENT  EFFICIENCY  PQ  FOR  RECTANGLE  AND  CORRESPONDING 
DATA  FOR  RIBBON  FROM  REFERENCE  11  (LANGMUIR) 


0.25 

0 

0.123 

0.122 

0 

0 

100 

.060 

.084 

0 

1,000 

.036 

.062 

0 

0 

10,000 

.011 

.035 

0 

0 

0.333 

0 

0.262 

0.244 

0.052 

0.098 

100 

.161 

.168 

.034 

.076 

1,000 

.113 

.125 

.017 

.063 

10,000 

.049 

.075 

.003 

.046 

0.50 

0 

0.436 

0.450 

0.287 

0.235 

100 

.316 

.350 

.231 

.186 

1,000 

.244 

.260 

.161 

.154 

10,000 

.134 

.157 

.066 

.114 

1 

0 

0.651 

0.666 

0.554 

0.470 

100 

.535 

.595 

.484 

.380 

1,000 

.452 

.503 

.384 

.320 

10,000 

.305 

.382 

.228 

.244 

2 

0 

0.789 

0.808 

0.722 

0.667 

100 

.691 

.734 

.650 

.574 

1,000 

.610 

.666 

.548 

.492 

10,000 

.466 

.565 

.388 

.389 

4 

0 

0.878 

0.894 

0.835 

0.800 

100 

.802 

.836 

.769 

.731 

1,000 

.733 

.783 

.676 

.660 

10,000 

.610 

.692 

.539 

.544 

10 

0 

0.953 

0.954 

0.928 

0.909 

100 

.903 

.918 

.874 

.864 

1,000 

.855 

.883 

.805 

.823 

• 

10,000 

.773 

.819 

.709 

.737 

20 

0 

0.982 

0.976 

0.962 

0.951 

100 

.948 

.954 

.923 

.921 

1,000 

.915 

.929 

.872 

.890 

10,000 

.862 

.882 

.803 

.834 

100 

0 

0.997 

0.995 

0.992 

0.989 

. 

100 

.982 

.988 

.970 

.978 

1,000 

.968 

.980 

.946 

.965 

10,000 

.948 

.960 

.912 

.942 

Figure  1.  -  Comparison  of  two-dimensional  flow  fields. 
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Local  Impingement  efficiency  at  stagnation  line. 
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Figure  4.  -  Local  impingement  efficiency  at  stagnation  line  calculated  from  equation  (20). 
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Figure  5.  -  Collection  efficiency  calculated  from  equation  (21). 
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